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Abstract 

We present a simple gluon production picture which is based on the McLerran-Venugopalan 
model and gluon BFKL evolution in relativistic heavy ion collision. Results for the multiplicity 
and transverse energy distribution in both the central and forward rapidity regions for gluon 
production in early stages of heavy ion collisions at the LHC are given. Finally, we provide a 
general qualitative discussion of the consequences of the forward rapidity behavior of produced 
gluons. 
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I. INTRODUCTION 



The description of the early stages of heavy ion collisions before hadronization is currently 
an interesting and evolving topic in heavy ion physics, and it is still not clear how good 
a description can be obtained from Quantum Chromodynamics(QCD). However, a better 
understanding of the properties of the initially produced gluons would shed some light on 
the early stages of the physics in heavy ion collisions. 

The McLerran-Venugopalan model [2] describes nuclear collisions at the Relativistic 
Relativistic Heavy Ion Collider (RHIC) in terms of the gluons in semi-classical non- 
abelian Weizsacker- Williams field generated by the valence quarks of relativistic nuclei. 
At RHIC energies, it appears to be a good approximation since the actual gluon occu- 
pation number(gluon phase space density) is large enough to be described by the classical 
non-abelian Weizsacker- Williams field. Specifically, we write the occupation number as 
fa = 2(%-i) dSfc = 2{lf-i) dy Sk x - At LHC energies, quantum evolution in the longitu- 
dinal momentum (BFKL evolution[l]) should be more important and increase the number 
of gluons by increasing the saturation momentum while the occupation number f g does not 
change muchj3j (we can see this either from the McLerran-Venugopalan model or from the 
discussion which follows.) 

The Balisky, Fadin, Kuraev and Lipatov (BFKL)[1] equation is an equation which de- 
scribes QCD hard scattering processes with a single hard transverse momentum scale and 
sums all terms like (a ln(l/x)) n (or (aY) n ). Explicit evaluation of this equation tells us 
that how nucleus evolves when it is boosted to a larger rapidity, and this gives the rapid- 
ity evolution of the nuclear gluon distribution. Moreover, it has been shown that the 
McLerran-Venugopalan model and fixed coupling BFKL evolution in the vicinity of the sat- 
uration boundary predicts geometrical scaling as shown below in Eq. (JSJ). We suppose that 
this geometrical scaling is the right description for the nuclear gluon distribution beyond 
the saturation momentum, that is when kj_ > Q s , at LHC energies. 

Based on a parton —on .nooeflflQ and toe oo,o r gl ass condensate^, as wen 



as a parton-hadron duality hypothesis, Kharzeev etal |9| have successfully made consistent 
calculation of hadron multiplicity in heavy ion collision, by assuming the produced gluons 
in the initial state hadronize into a corresponding number of hadrons in the final state. 
It appears that saturation and the Color Glass Condensate provide a reasonable physical 
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pictures for relativistic heavy ion collisions. Therefore, these tools have begun to be used 
to describe and understand the thermalization stages of the produced gluons right after 
their formation during the heavy ion collision. Thus, a description of the produced gluon 
multiplicity and transverse energy distribution throughout all rapidity regions seems to be 
more and more important, and this is part of the motivation of this work. 

In this paper, we are interested in the initial production processes of gluons in a heavy ion 
collision, which come from the freeing of the virtual gluons in colliding nuclei. The essential 
feature of this approach is that we choose our reference frame so that the produced gluon on 
which we are focusing has the rapidity y. Thus each of these two colliding nuclei have BFKL 
evolution in this particular frame. Based on the BFKL evolved gluon distribution in the 
nuclei and the approximate criteria J, ^| developed by Blaizot and Mueller, we calculate 
the distribution of freed gluons from four different possible regions which are shown in 
Eqs. (|Tl1 J) .(fTT )) .([T2 )) and (JT3J). In the end, we sum these contribution and obtain the total 
multiplicity and transverse energy distribution in Eqs. fl"5|) and (jl7)) for all produced gluons. 
In Eq. (|15|) . we find that the first and the third terms are similar to the results in ref. [^] 
while the second term is different since we use geometrical scaling in the region beyond 
the saturation momentum. We also calculate the average transverse momentum per gluon 
in forward rapidity region. The calculations show that the average transverse momentum 
increases a lot while the number of produced gluons decreases rapidly in the forward rapidity 
region which may lead to some interesting physics in that region. As we will discuss in detail 
later, we consider the consequence of this result for jet quenching, gluon thermalization and 
elliptic flow. One of the conclusions is that we might expect the jet quenching effect to be 
less significant in forward rapidity region. 

In section 2, we give a brief review of the McLerran-Venugopalan model, and derive the 
gluon distribution function in nuclei by following the interpretation of this model provided 
in ref . jlfll] . 

In section 3, we boost both of these two colliding nuclei (the projectile and the target) 
away from central rapidity region, and consider BFKL evolution |l| for both nuclei. 

In section 4, based on the framework established in above sections and the approximate 
criteria^, [ll | for gluon freeing, we calculate the rapidity dependence of the multiplicity and 
of the transverse energy of freed gluons. 

Finally, in section 5, some qualitative discussions regarding the consequences of the for- 
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ward rapidity behavior of produced gluons as well as a brief summary are provided. 



II. MCLERRAN-VENUGOPALAN MODEL 

Mueller and Kovechegov have shown that in the McLerran-Venugopalan model, the 
gluon distribution of the relativistic nucleus can be written as following: 

dxGjx, fcj) = N 2 - 1 r fx± / _ ^ Q 2 /4 x ik±.^ . , 

d 2 bd 2 k ± A^aNj ii 2 V ) ' 1 > 

in which the k± is the transverse momentum of the gluon, and the x — -p-p is the longitudinal 

momentum fraction of the gluon. The saturation momentum, Q s , is given by 

Q 2 S = ^^VW^ P xg(xA/x% (2) 

c 

in which the xg(x, 4/x 2 ) is just the gluon distribution of a single nucleon inside this nucleus. 
In the region k\ <C Q 2 , one can neglect the x 2 dependence of the Q 2 and evaluate the above 
gluon distribution integral 

dxGjx, kj) N 2 -l r dt y Ql N 2 c -l (Q 2 



d 2 bd 2 k ± 4n 3 aN r Ji t 4n 3 aN r \k 2 '' 




Moreover, for a large nucleus using BFKL evolution, it has been found |l2j, |13|, |14] that: 



dxGjx, fcj) _ N 2 - 1 1 - Ap (QT\ 
d 2 bd 2 k ± 4n 3 aN c 2 X (\o) n \k 2 ± ) ' U 

where A is defined by x'(^o) — — fi^> an d x(^o) = V'(l) ~~ IV^o) — 1^(1 ~~ -^o) with the 
tp = dln J^ defined as the logarithmic derivative of the gamma-function. It is likely that 
the factors 4 ^g ~^ and In (§r^J i n the gluon distribution formula can be taken as general 
results except for an overall factor which could be fixed by a comparison to RHIC data. 
The essential physical picture is that the gluon density per unit area and per unit of two- 
dimension transverse momentum, that is per unit of transverse phase space, is limited by 
(or saturates at) the product of a constant and logarithmic factor In (j^J with a upper 
momentum limit Q 2 . Therefore, one can write the gluon phase space density as 

dxG(x,k 2 ± ) kN 2 -l ln (Ql\ (5) 



d 2 bd?k\ ' c4n 3 aN r \ k 2 

2x(Ao 
1-Ao 



in which k is a 0(1) constant, and c = ^zt^ = 4.88. In practice, one can change above 



gluon distribution into a uniform distribution as -ga^j^ = ~ ^aN c ' s ^ nce ^ ne logarithmic 
factor disappears after integration over k±. 
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III. BFKL EVOLUTION AND GEOMETRICAL SCALING 



However, the simple McLerran-Venugopalan (MV) gluon distribution does not contain 
any evolution in Bjorken x which is necessary to explore the transverse energy spectrum. 
In order to extend our discussion from the central rapidity region to the forward rapidity 
region, we need to consider BFKL evolution of the saturation momentum. The leading order 
BFKL evolution for the saturation momentum is 



Q 2 S (A,Y)=Q 2 (A) 



exp 



2aN c xiXgl y 
it 1-Aq 



Numerical calculation shows that Aq 



(aY) 2 * 1 -^) 
0.372 and c 



2x(Ao) 
1-Ao 



(6) 



4.88. This leads to the 



saturation momentum Q 2 (Y) goes like exp [X S Y] with A» ~ 1 if taking a = 1/3. This value 
is much larger than the phenomenological approach es | l7l| which find that \ s is about 0.3. 
Fortunately, the next to leading order calculation |l8[ shows that A s effectively behaves 
like a constant around 0.30 since the next to leading correction slows down the evolution. 
Therefore, we shall take 

Q s (A,y) 2 = Q(A) 2 exp[0.30y] (7) 

in our numerical calculations. 

Nevertheless, the above simple MV gluon distributions Eqs. ©,© and (0) are only 
valid when k± << Q s , and they do not tell us how it behaves when the gluon transverse 
momentum k i is larger than the saturation momentum. Therefore, we need to consider the 



small-x evolution 



of the saturation momentum Q s and the gluon distribution in the 



vicinity of the saturation boundary. This gives a g eometrical scaling of the gluon distribution 



„,-,„„ „ , .0 L ™ i,i u™, 



16] 



saturation momentum, when k± > Q s , the evolved gluon distribution behaves like 



d?bd 2 k i 



k 



N 2 



ii (Ql 



1-Ao 



A-n: 3 aN r 



In 



(8) 



We have chosen the same overall constant to agree with Eq. ((SJ). 



IV. PRODUCED GLUON MULTIPLICITY AND TRANSVERSE ENERGY 



In this section we use the McLerran-Venugopalan model together with BFKL evolution 
to investigate behavior of produced gluons in all rapidity regions. Let us concentrate on 
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FIG. 1: Gluon distribution for two BFKL evolved nuclei, in which we denote the projectile(left- 
moving) nucleus saturation momentum as Q s (y + and target (right-moving) nucleus saturation 
momentum as Q s (—y + y). 

a simple case first, in which we consider the head on collision between the projectile (left- 
moving) nucleus and target (right-moving) nucleus where both of the nuclei have the same 
number of nucleon; therefore the physics is symmetric under y -<=>- — y. We also do not 
consider the impact parameter or centrality dependence for the moment. Thus for now we 
only consider the y > region. 

The essential feature of this approach is that we choose a reference frame in which the 
produced gluons has zero rapidity. The two colliding nuclei then are assumed to have 
saturation momenta in the BFKL evolution region. Suppose the total rapidity difference 
between the target and projectile is Y, then the saturation momenta for the target and 
projectile nuclei are Q s (—y + ¥) and Q s {y + ^) m this reference frame, respectively. Based 
on the BFKL evolved gluon distribution in the nuclei, we can determine the distribution of 
the freed gluons in our reference frame by using an approximate criteria^, 11] developed by 
Blaizot and Mueller just as one usually does in the central rapidity region in the center of 
mass frame. The approximate criteria states that the gluons which receive enough transverse 
momentum to push them onto mass shell will be freed during the collision. As far as we 
know, the typical momentum of the gluons inside the nuclei are Q s (—y + and Q s (y + y) 
respectively. Correspondingly, the necessary typical transfer momenta in order to free gluons 
are Q s (—y + and Q s (y + -r) respectively for these two nuclei as well. In this sense, we 
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FIG. 2: Gluon BFKL evolution in the projectile and target nuclei as well as two possible freeing 
scattering process during the collisions. 

can separate gluons inside the projectile and the target into four different according to these 
two saturation momenta as shown in Fig.((TJ. 



1. Gluons inside the target with the typical saturation momentum Q s (—y + 



a. Gluons inside the target with a transverse momentum k t < Q s (—y + (Noting that 
Q s (—y + y) < Qs(y + y) when y > 0) are labeled as 1 in Fig.(JTJ). They scatter with the 
gluons in the projectile nucleus which have a typical momentum Q s (y + \) due to the gluon 
saturation and BFKL evolution which are illustrated in Fig.(j2J). During the scattering, the 
virtual gluons which are originally bounded inside the target nucleus are freed according to 
the approximate criteria that we have stated above. Therefore, those gluons in part 1 inside 
the target can be freed (or more properly, they are produced), and we calculate their total 
number per unit rapidity per unit area from the definition = J®^ y+ 2 ^ d 2 k± dx d f^^ 
by integrating Eq. (jU ) 

dNt kN 2 c -l / , Y.\ 2 

Then, how are these produced gluons distributed in momentum space? Originally, we assume 
these gluons are almost uniformly distributed from QCD energy scale to the saturation 
momentum Q s (—y + ¥) according to Eq. (0J) and its simplified form where in Eq. the 



logarithmic term is replaced by a constant. During the freeing process, these gluons exchange 
transverse momentum in the scattering with gluons from the projectile which have a typical 



momentum Q s {y+^)- This means that the momentum space distribution of these gluons has 
been enlarged from Q s (— y + \) to Q s (y+^). Therefore, we presume that these gluons would 
be uniformly distributed in the momentum space from very small transverse momentum 
(could be close to zero) to the typical momentum of the projectile nucleus Q s (y + -j) scale 
which is the typical momentum the freed gluons can acquire after the freeing. Then one 
finally obtain: 

dN, _ kN 2 -l (Q s (-y + \)\ 2 



dyd 2 bd 2 k ± cA7r 3 aN c \ Q s (y + f ) J ' 
with k t < Q s (y + ~). Note that the distribution of the freed gluons has been broadened 
from kt < Q s (—y + to k t < Q s (y + ^) due to the scattering with the gluons from the 
projectile. 

b. Gluons inside the target with a transverse momentum k t obeying Q s (—y + ^) < k t < 
Qs(y + if) are labeled as 2 in Fig.(JTJ). They scatter with the gluons in the projectile nucleus 
which have a typical momentum Q s (y + \)- In this region, we believe that the distribution of 
the gluons inside the target obeys the geometrical scaling which we have shown in Eq. (JBJ). 
Therefore, these gluons can be freed and we assume they are uniformly distributed from 
very small transverse momentum to the typical momentum scale of the projectile nucleus 
Qs(y + y) scale after the freeing. The case is almost the same as the case a above except 
for the initial distribution of the gluons in the target. Then, we integrate over this region, 
and get the freed gluon distribution 



dNo k N 2 - 1 



dyd 2 bd 2 k ± c 4ir 3 aN ( 




:n) 



with kt < Q s (y + f-)- Actually, we neglect the logarithmic term again in Eq. (jSJ) be- 
fore arriving at the result in {} in above equation while the exact result would be 



\ \\„ ( Qs(v+t) 



.(-JH-t)/ W«(-9+f) 



In .\\ -1+1 



This simplification does not affect the 



physics much. In addition, we want to comment on the result that this is the largest 
contribution in the moderately large rapidity region due to the term which 

differs from the result in [l| although we both agree that this region gives largest contri- 
bution. This difference comes from the geometrical scaling, and we expect that heavy ion 
collisions at LHC might observe this scaling while at RHIC geometrical scaling may not yet 
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set in since the saturation momentum may not be large enough. 

c.Gluons inside the target with k t > Q s (y + y). In this region, the gluons are very few 
and the parton densities in the target are in the DGLAP evolution region. We believe this 
region does not contribute much. 

2. Gluons inside the projectile with the typical saturation momentum Q s (y + y). 

a. Gluons inside the projectile with a transverse momentum k t < Q s (—y+ j), labeled as 
3 in Fig, (fl)) . scatter with the gluons in the target nucleus which have a typical momentum 
Qs{~ V+ 2")- Therefore, these gluons can be freed and we assume they uniformly distributed 
from very small transverse momentum to the typical momentum scale of the target nucleus 
Qs(—y + j) after the freeing. Then one finds 

dN 3 _*AJ-1 



dyd 2 bd 2 k ± cA^aNc 

with k t < Q s (—y + y). Note that the range of the momentum remains the same because the 
typical momentum that the freed gluons scatter with is Q s (—y + y)- If we integrate over 
two-dimension phase space in Eqs. (|10p and (|12jh we get the same results with Kharzeev 
etal fl. 

b. Gluons inside the projectile with a transverse momentum Q s (—y + y) < k t < Q s (y + 
y) are labeled as 4 in Fig. (J!}. They scatter with the gluons in the target nucleus whose 
distribution is given by geometrical scaling. The problem here is very similar to that of case 
l.b. and we presume that the contribution for this part can be written as 

dN 4 _ g dN 2 
dyd 2 bd 2 k± dyd 2 bd 2 k± 

with < 5 < 1 and k t < Q s (y + j). 

c. Gluons inside the projectile with k t > Q s (y + y)- In this region the gluons are very few 
and the parton densities in the projectile are in the DGLAP evolution region. We believe 
that this region does not contribute much. 

Now, when we generalize our discussion to the very forward rapidity region, the factor x = 
^"^VS + 2 " ex P \.y] ^ s no l° n g er small. In this case we need to include a limiting frag mentation 
factor (1 — x) 4 which comes from the phenomenological quark counting rules 19J. 

Then, we add all contribution together and integrate over the two-dimension phase space 
to get the total multiplicity with the fragmentation factor taken into account as a large x 
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correction 



dN 
dy 



dyd 2 bd 2 k\ 



(1-x) 



TiR 2 ) {kQ 2 s (A, 




(14) 



-y + 77 



Y 
2 



1 



1 



1 



Qs(y + i) cy 



(15) 



In the above equation we have made the approximation that x = ^= exp [y] for the gluons 
with transverse momentum in the integral is so small that (1 — a;) 4 is always around 1 
until k± approaches the upper limits Q s (±y + ¥)■ 

Beginning with the approximate criteria of freeing gluon, we finally arrive at almost the 
same result as Kharzeev etal |2| except for the second term in above equation since we 
are using the geometrical scaling gluon distribution when k± > Q s .(We believe that the 
geometrical scaling distribution might be the right pict ure to describe gluon distribution 
beyond the saturation momentum at LHC energies [20].) Moreover, we also notice that 
the second term, which is the contribution from Eq. (fTTJl . is the largest term in forward 

Then we can go one 



rapidity region since it is enhanced by a factor of ' ^ s ^ y+ 2 - 1 



2A 



more step further, averaging over different transverse momenta to get the transverse energy 
distribution, 

l 



dE 7 
dy 



]T J d 2 bd 2 kj 



dNi 



i=l 

(k N? 



1 



4ir 3 aN r 



Q s {-y + 



+Qs(y + j) 



dyd 2 bd 2 k 

TlR 



xfk A 



(16) 



i - 



Q, 



1 + 5 

An 



Qs(y + j) 

2 • 



2A 



,Q s (~y 



i+i 



(17) 



In addition, we should fix the parameters which appear in the multiplicity and transverse 



energy distribution before we plot ^ and as a function of rapidity y in fig. 



dy 



We 



take the parameters to be a = 1/3 and N c = 3 according to convention, and we allow the 
parameter 5 to vary from to 1 in the plots. 
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FIG. 3: Predictions of the multiplicity transverse energy distribution and average trans- 
verse momentum per gluon (k±(y)} for the V^ = 5500GeV lead-lead head-on collision at LHC 
when using Q s = l.95GeV and a = 1/3, respectively. 



Noting that exp[^4p] = where M nuc is the rest mass of one of the colliding 

nuclei, we can write the Eq.( 



as 



Q 2 (A,y+-) = Q 2 So (A) 
Then, phenomenologically we take: 

R = r A^, 
Q 2 S {A) = (l.lGeV) 2 



exp [X a y] . 



(A 



A 
V197 



(18) 

(19) 
(20) 



\130GeV J 

where = 1.2 fm is obtained from the empirical nuclear radius formula and Q s 1,4=197 = 
l.lGeV 3] is an estimation based on RHIC Au-Au collision at ^/sq = 130GeV by using 
Eq. (J2J) and taking p = i^ R3 = t^3, xg(x, Q 2 ) = 3.3. Therefore, we predict that the satura- 



3 nT o 



tion momentum at LHC is around 1.95GeV at ^/s = 5500GeV in the central rapidity region 
and ttR 2 = 158 fm 2 for lead nuclei from Eqs. (fl"9|) and 1)20)1 . Noting that the geometrical 
scaling contribution in the multiplicity Eq. (|15j) vanishes in the y = central rapidity region 
means that \ y= o should match RHIC data. Thus, we can fix the ratio of parameters 
k and c by using the multiplicity data in the central rapidity region from PHOBOS Q 



which provides 



dy 



^/s=130GeV 

y=o 



547 ± 55 or 



dN 
dy 



,/i=130GeV 
y=0 



821 ± 83 (Here we assume 



dN I 

dy \y=° 



3 dN ch 
2 dy 



j /=0 in which | y=0 is the multiplicity of charged particles, and we also 
neglect the difference between the rapidity y and pseudorapidity 77.). This helps us to fix 
I = 2"jj4 which agrees with Mueller's conjecture 111] ( Mueller suggests that the total number 
of the produced gluons is proportional to the total number of the gluons in the initial nuclear 
wave-function and the proportionality coefficient e is called gluon liberation coefficient.) and 

25j that the gluon liberation coefficient e = — = ttIt (this re- 



previous studies [2z 



23, 



24, 



dN 



lation can be extracted from Eq. (JT5j) via the definition d2My \ y=0 
rapidity region.) is an order 1 coefficient. 

Then with fixed coefficients and Eq. (|15l) . we are able to predict that: 



- N 2 c ~l 
'4n 3 aN c 



ttQ 2 at central 



dNch 
dy 



v/i=200GeV 
y=0 



630 which agrees with PHOBOS data[26] 



dN ch 

dy 



v^=200GeV 
y=0 



650 ± 35 (syst) and 



dNch 
dy 



^=5500GeV 
y=0 



1780 for the lead- lead collision in the LHC energies, which is close to 



Kharzeev's prediction 



19] 



dN ch 
dy 



N /s=5500GeV 

y=o 



1750 - 2100). 



Finally, for each rapidity, we calculate the average transverse momentum per gluon as 

dE T 

(k±(y)) = -Sfr (21) 

dy 
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and plot (k±(y)) also in Fig. Q. 

We find that the average transverse momentum per parton in our calculation (as well as 
that in Kharzeev's paper if one generalize it and calculate the average transverse momentum) 
increases a lot as we move to the forward rapidity region. (When y > 6, (k±(y)) starts to 
decrease because of the nonlinear effect of the limiting fragmentation factor (1 — x) 4 . We do 
not expect our discussion to hold in this region since x again is no longer small when y > 6, 
even at LHC energies.) 



However, what we find is very different from the forward rapidity result [27( from RHIC 
for the final state hadrons where the average transverse momentum per hadron is almost 
constant as rapidity increases. There might be two reasons to explain this discrepancy. First, 
at RHIC the typical Bjorken x is less than 10~ 2 , and becomes comparable to 1 in the forward 
rapidity region. Also, the saturation momentum for one of the colliding nuclei becomes very 
small. This makes the saturation picture not applicable to the RHIC data in that region. The 
second reason is that there exist several intermediate processes (e.g., thermalization.2j| and 
hadronization) between the initial produced gluons and final measured particles. This may 
indicate that the direct use of the parton-hadron duality hypothesis might be questionable 
and dangerous in this case. In order to get the right picture from initially produced partons 
to the final state particles, we need to consider the intermediate processes like thermalization 
which distorts the transverse momentum distribution of produced partons. 

Finally, we sketch how to deal with the case where the projectile and target have different 
number of nucleons(e.g. D-Au collision). As one can see, ^ reaches its maximum when two 
saturation momentum are the same. This happens at y — in the A\ = A2 case. While in 
the D-Au collision, the critical rapidity can be obtained by setting 

Qs(A u -y + ^) 2 = Q,{Ai? M= \ S exp[-A s2/ ] 



Q S (A 2 , V+\f = Q*(M) 2 S exp [X s y] (22) 



and noting that the saturation momentum is proportional to A 



3 . 
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V. CONCLUSION 



In the following, we would like to qualitatively discuss a little bit about the consequences 
of the above results in forward rapidity region. As we have seen (presuming 5 is small 
compared to 1), the number of the produced gluon and the transverse energy in per unit 
rapidity decrease considerably in forward rapidity region, while the transverse energy per 
gluon per unit rapidity increases a lot in forward rapidity region. This leads to some changes 
compared to the central rapidity region. 

1. Jet quenching: Jet quenching is a very important signal for the existence of the 
quark gluon plasma (QGP) in relativistic heavy ion collisions. After the hard scattering 
between two partons, each of the two scattered partons evolves into a jet which consists of 
a leading particle together with a number of associate hadrons around the leading particle 
during hadronization. However, in the presence of dense strongly interacting matter (QGP) 
surrounding the hard scattering site, the scattered partons have induced gluon radiation 
due to multiple scattering with the plasma and this causes the jet to lose energy. This 
energy loss causes a suppression of the observed leading hadron spectrum. This is the jet 
quenching. However, since the number of the produced gluons is much less in the forward 
rapidity region than in the central rapidity region, the induced energy loss is smaller since 
bremsstrahlung gluon radiation and multiple scattering are expected to be less important. 
Therefore, we might expect that the jet quenching effect would be less significant in the 
forward rapidity region. However, in the very large rapidity region our discussion breaks 
down because of two reasons: first, the saturation is no longer applicable to the region where 
x is large; the other reason is that most of the valence quarks are in this region after the 
collision. The valence quarks may give a significant contribution to partonic density and 
increase the induced strong interaction which then might result in more energy loss. 

2. Gluon thermalization: Baier eta/.j^ have discussed thermalization of the partons, 
after they are produced and before hadronization, in the bottom-up thermalization scenario. 
In the above picture, thermalization time can be expected to be considerably longer in the 
regions away from the central rapidity as compared to the thermalization a t y = 0. 



29 



30 



31 



is intro- 



3. The elliptic flow: The concept of the elliptic flow parameter t>2 
duced to reflect the anisotropy of the produced gluons that due to the spatial geometrical 
azimuthal asymmetry in the non-central collision. With the similar picture and discussion 
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for the jet quenching, we can see that the elliptic flow is usually generated at the partonic 
level prior to the hadronization. At RHIC, data exhibits a very strong elliptic flow and 
saturate the hydrodynamic limit. This suggests that very rapid and almost complete local 
kinetic equilibrium is reached. At LHC, elliptic flow is again expected to saturate the hy- 
drodynamic limit. In the forward rapidity region at LHC energies due to the decrease of 
produced gluons, we may expect a smaller elliptic flow. 

In the saturation picture, we apply the McLerran-Venugopalan model and the geomet- 
rical scaling to the computation of the multiplicity and transverse energy distribution for 
heavy-ion collisions at LHC. Results for the multiplicity and transverse energy distribution 
in the central rapidity region agree well with the RHIC experimental data and previous 
studies. Then we make predictions for the results of the produced gluons at the early stage 
of heavy ion collisions in both central and forward rapidity regions at the LHC energy scale. 
Noticing that hadronization and thermalization can distort both the multiplicity and trans- 
verse energy distribution, we need to consider these effects before we compare the results of 
the produced gluons with the data of the final state particles. 
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